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The method involves evaluating sound from the 
target received directly and via reflections. The 
hardware includes a submerged receiver with a 
number of electroacoustic transducers. 
The electrical output signals off the transducers 
are added directionally selectively to form group 
signals and transition time differences between 
direct and indirect paths determined and 
compared with modeled differences with various 
target position assumptions until agreement is 
achieved. 
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(54) Verfahren zum Bestimmen der Zielposition eines schallabstrahlenden Ziels 



(57) Verfahren zum Bestimmen der Zielposition ei- 
nes schallabstrahlenden Ziels unter Auswertung des 
vom Ziel abgestrahlten Schalls, der auf direktem Aus- 
breitungsweg und auf durch Reflexlonen entstehenden 
indirekten Ausbreitungswegen von einem Schallemp- 
fanger mtt einer Mehrzahl von elektroakustischen 
Wandlern empfangen wird. Zu einer verbesserten Ziel- 
positionsbestimmung von breitbandig Schall abstrah- 
lenden Zielen werden richtungsselektive Gruppensi- 
gnale gebildet und Laufzeitdifferenzen zwischen dem 
direkten Ausbreitungsweg und indirekten Ausbreitungs- 
wegen gemessen. Mittels eines Modells der Schallaus- 
breitung fur angenommene Zielpositionen werden zu- 
gehorige, gleiche Modell-Laufzeitdifferenzen bestimmt 
und mit den MeB-Laufzeitdifferenzen verglichen. Die 
Bestimmung der Modell-Laufzeitdifferenzen mitjeweils 
geanderter Zielpositionsannahme wird solange wieder- 
holt, bis der Vergleich einen vorgegebenen Grad der 
Ubereinstimmung ergibt. Die zum Erreichen dieses 
Ubereinstimmungsgrads getroffene Zielpositionsan- 
nahme wird als echte Zielposition ausgegeben. Fehler 
durch den Dopplereffekt werden beim Modellieren von 
M ustersignalen fur die Bestimmung der MeB-Laufzeitdif- 
ferenzen beseitigt und eine Geschwindigkeitsbestim- 
mung ermoglicht. 
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Beschreibung 

[0001 ] Die Erfindung betrifft ein Verfahren zum Bestimmen der Zielposition eines Schall ins Wasser abstrahlenden 
Ziels der im Oberbegriff des Anspruchs 1 definierten Gattung. 

s [0002] Bei einem bekannten Verfahren dieser Art (DE 37 05 695 C1 ) wird der mit kleinen elektroakustischen Wand- 
lern, dort Sensoren genannt, ausgestattete Schallempfanger, dort Sonarempfanger genannt, in einer ersten Wasser- 
tiefe angeordnet und ein erster Naherungswert der Laufdifferenzen des Schalls zwischen dem direkten Ausbreitungs- 
weg und einem indirekten Ausbreitungsweg mit Oberflachenreflexion zu einer ersten Zeit bestimmt. AnschlieBend 
bewegt sich der Schallempfanger mit einer Geschwindigkeit in eine zweite Wassertiefe und ein zweiter Naherungswert 

10 der Laufzeitdifferenz des Schalls zwischen dem direkten Ausbreitungsweg und einem indirekten Ausbreitungsweg mit 
Oberflachenreflexion zu einer zweiten Zeit bestimmt. Auf der Grundlage der beiden Naherungswerte fur die Laufzeit- 
differenzen, einem Naherungswert fur die Geschwindigkeit des Schallempf angers und dem Zeitintervall zwischen er- 
ster und zweiter Zeit wird ein Naherungswert fur die Zielentfemung geschatzt und daraus ein Naherungswert fur die 
Zieltiefe abgeleitet. 

is [0003] Bei einem unter dem Begriff "Pingstealing" bekannten Verfahren, das ein aktiv Schallimpulse abstrahlendes 
Ziel - wie es z.B. ein aktiv ortendes Torpedosonar eines Torpedos darstellt - voraussetzt (US 4 312 053), wird mittels 
eines getauchten Schallempfangers der Einfallswinkel der akustischen Energie innerhalb eines Schallimpulses als 
Funktion der Zeit gemessen. Diese Daten werden dann in Mehrwege-Einfallswinkel und Mehrwege-Laufzeitdifferenzen 
umgesetzt und damit ein mathematisches Gleichungssystem zur Berechnung der Zielentfemung und Zieltiefe gelost. 

20 [0004] Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren der eingangs genannten Art zu schaffen, das auch 
mit vom Ziel breitbandig abgestrahlter akustischer Energie, wie z.B. Antriebsgerausche eines sich bewegenden Ziels, 
ohne zusatzliche Bedingungen, wie z.B. eine Horizontal- und Vertikalbewegung des Schallempfangers, insbesondere 
in Flachwassergebieten sehr genaue Ergebnisse bei der Zielvermessung liefert. 
[0005] Die Aufgabe ist erfindungsgemaB durch die Merkmale im Anspruch 1 gelost. 

25 [0006] Das erfindungsgemSBe Verfahren hat den Vorteil, dal3 durch die Bildung der Gruppensignale zum Zwecke 
eines richtungsabhangigen Anhebens des Nutz-Stor-Verhaltnisses (SNR) und durch die Auswertung der in den Grup- 
pensignalen enthaltenen Laufzeitdifferenzen zwischen unterschiedlichen Ausbreitungswegen mittels Vergleich mit ei- 
ner Modellannahme der Schallausbreitung auch bei geringen Wassertiefen die Zielposition aufgrund der vom Ziel 
abgestrahlten Gerausche zuverlassig genau bestimmt wird. Bereits mit einfachen geometrischen Ausbreitungsmodel- 

30 len fur den Schall kann die Zielposition mit sehr hoher Wahrscheinlichkeit bestimmt werden. Mit der Kenntnis der 
Zieltiefe laBt sich dariiber hinaus der azimutale Peilwinkel des Ziels mit einem berechenbaren Elevationswinkel des 
Schalleinfalls korrigieren und so das Toleranzgebiet der Zielposition einengen. Das Verfahren ist auch in der Lage, ein 
Ziel zu detektieren und die Zielpeilung zu bestimmen, wenn die richtungsselektive Bildung der Gruppensignale konti- 
nuierlich uber den gesamten Azimut durchgefiihrt wird und auf alle Gruppensignale die erfindungsgemaBe Signalver- 

35 arbeitung angewandt wird. Selbstverstandlich ist das Verfahren auch bei Zielen anwendbar, die aktiv Schallimpulse 
senden. 

[0007] ZweckmaRige Ausfuhrungsformen des erfindungsgemaBen Verfahrens mit vorteilhaften Weiterbildungen und 
Ausgestaltungen der Erfindung ergeben sich aus den weiteren AnsprOchen. 

[0008] Der Abgleich zwischen Modell und Messung kann auf verschiedenen Ebenen stattfinden, wobei der Grad der 
40 Obereinstimmung ein Maf3 fur die Aufenthaltswahrscheinlichkeit des Ziels in der angenommenen Zielposition ist, die 
durch Peilung, Entfernung und Zieltiefe kennzeichnend ist. 

[0009] Bei einer bevorzugten Ausfuhrungsform der Erfindung wird ein direkter Vergleich der Modell- und MeB-Lauf- 
zeitdifferenzen durchgefiihrt, indem als Kriterium die Summe der minimalen Abstandsquadrate der Modell-Laufzeit- 
differenzen zu ihren nachsten Nachbam in der Menge der MeR-Laufzeitdifferenzen verwendet wird. Dervorgegebene 
■*5 Grad der Obereinstimmung ist dann erreicht, wenn die Summe dieser minimalen Abstandsquadrate innerhalb eines 
Konfidenzintervalls liegt. 

[0010] GemaB einer vorteilhaften Ausfuhrungsform der Erfindung werden dabei die minimalen Abstandsquadrate 
vor Summenbildung mit einem Wichtungsfaktor Wj multipliziert, wobei als Wichtungsfaktor W; eine Potenz der Anzahl 
n B ,n D mit den Reflexionsfaktoren R Q im Ausbreitungsweg i gemaB 

50 

w, = R B B -Ro° 0) 

benutzt wird. Der Reflexionsfaktor R B am Boden B und der Reflexionsfaktor an der Oberflache 0 eines durch Boden 
55 B und Oberflache O begrenzten Schallkanals wird plausibel zwischen 0 und 1 gewahlt. Durch diese Wichtung wirken 
sich stark gedampfte und nicht identifizierbare MeB-Laufzeitdifferenzen hoherer Ordnung nicht verfalschend auf die 
Aufenthaltswahrscheinlichkeit aus. 
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[001 1] Die Identifizierung der Laufzeitdifferenzen in jedem aus den elektrischen Ausgangssignalen von i elektrischen 
Wandlern gebildeten Gruppensignal, kann nach verschiedenenMethoden erfolgen. Die Gruppensignale, auch Beams 
genannt, werden durch koharente Uberlagerung derentsprechend der Einfallsrichtung des Schalls laufzeitverzogerten 
Ausgangssignale gebildet. 

[0012] GemaB einer vorteilhaften Ausfuhrungsform der Erfindung wird zur Laufzeitanalyse die Autokorrelationsfunk- 
tiondesGruppensignalsgemessen, und die MeB-Laufzeitdifferenzen fur die unterschiedlichen Ausbreitungswege wer- 
den zu vorgebbaren MeBzeitpunkten durch Auswertung der Maxima der Autokorrelationsfunktion gewonnen. Zwischen 
zwei MeBzeitpunkten liegt ein Zeitintervall, innerhalb dessen sich das Ziel weiterbewegt hat. Da die Autokorrelations- 
funktion eines breitbandigen Signals mit der Bandbreite B Nebenmaxima mit abnehmender Amplitude im Abstand 
At=1/B besitzt, lassen sich in der Autokorrelationsfunktion benachbarte Hauptmaxima, die die MeB-Laufzeitdifferenzen 
kennzeichnen, nurdann voneinandertrennen, wenn die Amplituden der Nebenmaxima ausreichend schnell abfallen, 
was bedingt, da3 das Signal geeignet breitbandig sein muB. 

[0013] GemaB einer bevorzugten Ausfuhrungsform der Erfindung wird die Identifizierung der im Gruppensignal ent- 
haltenen Laufzeitdifferenzen mittels einer Cepstralanalyse vorgenommen. Das reelle Cepstrum ist definiert als Fou- 
riertransformierte des Logarithmus des Absolutwerts des Autospektrums eines Signals 



wobei Xfe^t) die Fouriertransformierte des Signals darstellt. Zur Definition des Autospektrums vgl. [2], Seite 50 ff. Die 
im Autospektrum enthaltenen Laufzeitdifferenzen lassen sich in der Fouriertransformierten des Spektrums, also im 
Cepstrum, anhand der dort auftretenden Maxima identifizieren. Zur Identifizierung der Maxima spielen die Langen der 
Fouriertransfonnation fur das Autospektrum und dessen Fouriertransformierten eine entscheidende Rolle. Die Lange 
des Autospektrums soil daher moglichst groB sein, darf aber die zeitliche Koharenzlange des Signals nicht uberschrei- 
ten, da dann die Maxima im Cepstrum wieder verschwinden. 

[0014] GemaB einer alternativen Ausfuhrungsform der Erfindung kann bei der Identifizierung der im Gruppensignal 
enthaltenen Laufzeitdifferenzen durch Bestimmung des Cepstrums des Gruppensignals der Vergleich der Modell-Lauf- 
zeitdifferenzen mit den MeB-Laufzeitdifferenzen in der Weise vorgenommen werden, daB mit den Modell-Laufzeitdif- 
ferenzen ein Modell-Cepstrum gebildet wird und Modell-Cepstrumund gemessenes Cepstrum unmittelbarmiteinander 
verglichen werden. 

[0015] Als Modelle fur die Schallausbreitung zur Bestimmung der Modell-Laufzeitdifferenzen aufgrund von Zielpo- 
sitionsannahmen werden von der Vielzahl bekannter Modelle fur die Mehrwegeausbreitung des Schalls in einem 
Schallkanal die geeigneten Modelle herangezogen. Solche Modelle sind beispielsweise aus [1] ablesbar. 
[0016] GemaB einer vorteilhaften Ausfuhrungsform der Erfindung wird ein einf aches geometrisches Modell in einem 
Flachwassergebiet mit der Wassertiefe z w bei ebenem Gewasserboden und konstantem Schallgeschwindigkeitsprof il 
zugrundegelegt. Bei einer bekannten Tiefe Z E des Schallempfangers und einer angenommenen Zieltiefe z s und azi- 
mutalen Zielentfemung x wird die Schallaufzeit t^ fur k ausgewahle Ausbreitungswege gemaB 



berechnet, wobei c die konstante Schallgeschwindigkeit und z k der Weg ist, den der Schall in jedem Ausbreitungsweg 
k ausschlieBlich in vertikaler Richtung zurucklegt. 

[0017] Mit einem solchen einfachen geometrischen Modell werden mit geringem Rechenaufwand bereits ausrei- 
chend genaue Ergebnisse in der Vermessung der Zielposition eines detektierten Zieis erzielt. 
[001 8] GemaB einer vorteilhaften Ausfuhrungsform der Erfindung wird bei festgestellter Zielposition die radiale Ge- 
schwindigkeit des Zieis, d.h. die auf den Empfanger gerichtete Geschwindigkeitskomponente des Ziels, und die Fre- 
quenz des vom Ziel abgestrahlten Schalls ermittelt. Hierzu werden die bei der Bestimmung der Modell-Laufzeitdiffe- 
renzen notwendigerweise berechneten Schall-Laufzeiten fur den direkten Ausbreitungsweg des Schalls und fur die 
ausgewahlten indirekten Ausbreitungswege des Schalls zur festgestellten Zielpositionen und das Gruppensignal, in 
dem das Ziel detektiert worden ist, aufgerufen. Zu Zeitpunkten, die den aufgerufenen Schall-Laufzeiten entsprechen, 
werden aus dem Gruppensignal Frequenzspektren gebildet. Aus den Frequenzspektren werden Empfangsfrequenzen 
ausgelesen, die Differenzen der Empfangsfrequenzen gebildet und den aufgerufenen Schall-Laufzeiten zugeordnet. 
Fur eine angenommene, vom Ziel abgestrahlte Frequenz und fur eine angenommene, radiale Geschwindigkeit des 
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Ziels werden zu jeder der aufgerufenen Schallaufzeiten Dopplerfrequenzverschiebungen, kurz Dopplerverschiebun- 
gen, berechnet und die Differenzen der Dopplerverschiebungen gebildet. Die Empfangsfrequenz-Differenzen und die 
Dopplerverschiebungs-Differenzen werden miteinanderverglichen, und die Berechnung der Dopplerverschiebungen 
mit jeweils geanderter Zieifrequenz und/oder radialer Zielgeschwindigkeit und der Vergleich mit den Empfangsfre- 
quenz-Differenzen wird solange wiederholt, bis sich ein vorgegebener Grad der Ubereinstimmung ergibt. Die zum 
Erreichen dieses Obereinstimmungsgrads getroffenen Annahmen von Zieifrequenz und radialer Zielgeschwindigkeit 
werden als wahre Zieifrequenz und wahre radiale Zielgeschwindigkeit ausgegeben. 

[0019] Diese Weiterbildung des erfindungsgemaBen Verfahrens hat den Vorteil, da3 durch rein passive Detektion 
zusatzliche Zielparameter erhalten werden, die wertvollen AufschluB uber eine mogliche Bedrohung durch das Ziel 
geben. Mit der Bestimmung der Zieifrequenz kann das Ziel klassifiziert und mit der radialen Zielgeschwindigkeit die 
Geschwindigkeit festgestellt werden, mit der sich das Ziel auf den Empfanger zubewegt. Dies ermoglicht z.B. eine 
fruhzeitige Torpedodetektion und -wamung und die Einleitung von AbwehrmaBnahemen aufgrund der Kenntnis der 
Torpedoangriffsgeschwindigkeit. 1st das Aktivsonar des Torpedos aktiviert, so ergeben sich in den Frequenzspektren 
des Gruppensignals signifikante Frequenzlinien, die der dopplerverschobenen Sendefrequenz des Aktivsonas ent- 
sprechen, wobei die Dopplerverschiebungen uber die verschiedenen Schallausbreitungswege unterschiedlich sind. 
[0020] Bei nicht aktiviertem Torpedosonar wird der Torpedo aufgrund seines spezifischen Antriebsgerausches de- 
tektiert und seine Zielposition aus MeB-Laufzeitdifferenzen und Modell-Laufzeitdifferenzen bestimmt. In den Frequenz- 
spektren des zur Zielposition gehorigen Gruppensignals treten signifikante hochfrequente Frequenzbander auf, deren 
Mittenfrequenz mit ausreichender Genauigkeit Informationen iiber die darin enthaltene Dopplerverschiebung zur Be- 
stimmung der radialen Zielgeschwindigkeit geben. 

[0021] GemaB einer bevorzugten Ausfuhrungsform der Erfindung wird die Dopplerverschiebung im direkten Aus- 
breitungsweg des Schalls als Produkt aus angenommener Zieifrequenz und angenommener radialer Zielgeschwin- 
digkeit dividiert durch die Schallgeschwindigkeit im Wasser bestimmt. Die Dopplerverschiebungen fur die ausgewahl- 
ten indirekten Ausbreitungswege werden aus dem Verhaltnis der Laufzeit uber dem direkten Ausbreitungsweg zur 
Laufzeit uber den jeweiligen indirekten Ausbreitungsweg multipliziert mit der Dopplerverschiebung fur den direkten 
Ausbreitungsweg berechnet. 

[0022] GemaG einer weiteren vorteilhaften Ausfuhrungsform der Erfindung wird die Positionsbestimmung des Ziels 
verbessert und die Bestimmung der Geschwindigkeit des Ziels ohne eine Frequenzanalyse moglich. Der sich vom Ziel 
ausbreitende Schall wird auf dem direkten Ausbreitungsweg zum Empfanger am wenigsten gedampft, da sein Weg 
am kurzesten ist und Verluste durch Reflexionen nicht auftreten. Durch die Geschwindigkeit des Ziels wird der Schall 
abhangig vom Ausbreitungsweg bei einer Annaherung des Ziels zeitlich unterschiedlich entsprechend dem Doppler- 
effekt gestaucht, so daG im Abstand der Laufzeitdrfferenzen die Dopplerverschiebung des abgestrahlten Schalls un- 
terschiedlich und zum direkt empfangenen Schall am groBten ist. Diese gegenuberder Direktausbreitung unterschied- 
liche zeitliche Stauchung Oder Dehnung des Schalls am Empfangsort wird erfindungsgemaG bei der Laufzeitbestim- 
mung des Schalls fur jede Posrtionsanderung des Ziels berucksichtigt. 

[0023] Im Modell wird ein von der Geschwindigkeit abhangiger Modell-Stauchungsfaktor je Ausbreitungsweg fur die 
nachste angenommene Position des Ziels bereitgestellt. Zur Schitzung der nachsten angenommenen Zielposition 
wird zu der bereits bestimmten Zielposition ein wahrend eines Zeitintervalls mit einer angenommenen Geschwindigkeit 
zuruckgelegtes Wegintervall hinzugefugt und dazu die Modell-Laufzeitdifferenzen und ihre zeitlichen Ableitungen be- 
stimmt, die die Modell-Stauchungsfaktoren bilden. Sie entsprechen der relativen Dopplerverschiebung und werden fur 
unterschiedliche Geschwindigkeiten aus der vorherigen Position fur angenommene Positionen je Ausbreitungsweg 
aus der zeitlichen Anderung der Modell-Laufzeitdifferenzen berechnet und den Modell-Laufzeitdifferenzen zugeordnet. 
[0024] Zu dem MeBzeitpunkt, zu dem die Position des Ziel bereits ermittelt wurde, wird das Gruppensignal mit dem 
groBten Pegel innerhalb eines Zeitintervalls gespeichert und dem direkten Ausbreitungsweg zugeordnet. Seine Dopp- 
lerverschiebung bleibt konstant wenn wahrend der nachsten beiden Zeitintervalle die Geschwindigkeit des Ziel auf 
dem direkten Ausbreitungsweg als konstant angenommen wird. Zum Erstellen von Mustersignalen wird das abgespei- 
cherte Gruppensignal jeweils mit dem Modell-Stauchungsfaktoren moduliert, indem der zeitliche Verlauf durch unter- 
schiedliche Ausspeicherzeiten gestaucht oder gedehnt wird. Fur jede angenommene Geschwindigkeit und Zielposition 
wird je Ausbreitungsweg ein Mustersignal erstellt. 

[0025] Zum nachsten MeBzeitpunkt wird das momentan empfangenen Gruppensignal mit den Mustersignalen, die 
aus dem Gruppensignal des vorherigen MeBzeitpunkts gewonnen wurden, kreuzkorreliert. Jede Kreuzkorrelations- 
funktion, die im Zeitoder Frequenzbereich erfolgen kann, liefert bei dem mit dem Mustersignal Qbereinstimmenden 
Ausbreitungsweg eine MeG-Laufzeitdifferenz, die der Modell-Laufzeitdifferenz entspricht. Zugleich ist die angenom- 
mene Geschwindigkeit gleich der Geschwindigkeit des Ziels, mit der es die neue Position eingenommen hat, wenn 
samtliche MeG-Laufzeitdifferenzen und Modell-Laufzeitdifferenzen fur alle Ausbreitungswege iibereinanderliegen und 
der Vergleich zwischen Messung und Modell die kleinste Abweichung aufweist. 

[0026] Der Vorteil der erfindungsgemaBen Weiterbildung besteht darin, daG die auch durch die Ambiguityfunktion 
bekannten Fehler bei der Laufzeitbestimmung, die durch die Stauchung Oder Dehnung des empfangenen Gruppen- 
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signals gegenuber dem abgestrahlten Schall auftreten, durch die Bereitstellung der Mustersignaleverkleinert werden. 
Durch die zeitliche Dehnung oder Stauchung des abgespeicherten Zeitstucks des Gruppenslgnals, das dem dlrekten 
Ausbreitungsweg zuzuordnen ist, wird ein Mustersignal geschaffen, das auf den uber den Umweg empfangenen Schall 
frequenz- und zeitmaBig angepaBt ist. Nurbei Ubereinstimmung liefern die Kreuzkorrelationsfunktionen Maxima. Wer- 
den Impulse empfangen, wird durch die erfindungsgemaBe Erstellung der Mustersignale die Dehnung und Stauchung 
der Impulslange je Ausbreitungsweg bereits berucksichtigt. 

[0027] Der Vorteil des erfindungsgemaBen Verfahrens besteht darin, daB breitbandig ohne Frequenzmessung allein 
aus der mit dem Dopplereffekt erklarbaren zeitlichen Stauchung oder Dehnung der ubertragenen Schaliwelle die Ge- 
schwindigkeit des Ziels ermittelt werden kann, da die Modell-Stauchungsfaktoren fur die unterschiedlichen Ausbrei- 
tungswege abhangig vom Schallkanal und einer vorgegebenen Geschwindigkeit berechnet werden. 
[0028] Fur eine Klassifizierung des Ziels besteht der Vorteil, aus der relativen Dopplerverschiebung, die der Modell- 
Stauchungsfaktor angibt, und der Geschwindigkeit des Ziels die Frequenz oder das Frequenzband des vom Ziel ab- 
gestrahlten Schalls ohne Frequenzanalyse zu ermitteln. 

[0029] Eine weitere Verbesserung, insbesondere bei einer Annaherung des Ziels, besteht darin, daB ab einer vor- 
gebbaren Zielentfernung von einer breitbandige Auswertung des Gruppensignals auf eine Auswertung innerhalb eines 
Frequenzbreichs mit vorgebbarer Mittenfrequenz und Bandbreite umgeschaltet werden kann und dadurch das NutzV 
Storverhaltnis verbessert wird. Insbesondere im mittleren und oberen Frequenzbereich, der im Nahbereich trotz der 
hoheren Dampfung gut empfangen werden kann, ist die Dopplerverschiebung fur die unterschiedlichen Ausbreitungs- 
wege signif ikant unterschiedlich, so daB die Bestimmung der Maxima der Kreuzkorrelationsf unktion fiir die unterschied- 
lichen Ausbreitungswege eindeutige Ergebnisse liefert und damit die Positions- und Gerschwindigkeitsbestimmung 
verbessert wird. Weiterhin besteht der Vorteil, daB durch die Bandbreite die Auflosung der Entfernungsmessung und 
durch die Integrationszeit bei der Kreuzkorrelationsfunktion die Auflosung der Geschwindigkeitsmessung abhangig 
von einer Falschalarmrate eingestellt werden kann. Die Entfernungsauflosung ist umso besser, je groBerdie Bandbreite 
ist, wahrend die Geschwindigkeitsauflosung mit steigender Integrationszeit wachst. Zumindest wahrend der Integra- 
tionszeit und des Zeitintervalls muB die Geschwindigkeit konstant und das Qbertragungsverhalten des Schallkanals 
stabil bleiben. 

[0030] Die Erfindung ist anhand eines in der Zeichnung illustrierten Ausfuhrungsbeispiels im folgenden naher be- 
schrieben. Es zeigen: 

Fig. 1 ein Blockschaltbild einer Schaltungsanordnung zur Erlauterung des Verfahrens zum Bestimmen der Zielpo- 
sition eines schallabstrahlenden Ziels, 

Fig. 2 ein Diagramm einer Autokorrelationsfunktion eines Gruppensignals zur Bestimmung der MeB-Laufzeitdiffe- 
renzen, 

Fig. 3 ein Diagramm eines Cepstrums des gleichen Gruppensignals zur Bestimmung der MeB-Laufzeitdifferenzen, 

Fig. 4 eine grafische Darstellung eines geometrischen Schallausbreitungsmodells mit Berechnung der Schallauf- 
zeiten in den einzelnen Ausbreitungswegen, 

Fig. 5 ein Blockschaltbild einer gegenuber Fig. 1 modifizierten Schaltungsanordnung, 

Fig. 6 ein gegenuber dem Blockschaltbild gemaB Fig. 1 erweitertes Blockschaltbild zur Erlauterung des Verfahrens 
zum Bestimmen der Zielparameter: radiale Zielgeschwindigkeit und Zielfrequenz nach festgestellter Zielpo- 
sition und 

Fig. 7 ein gegenuber dem Blockschaltbild gemaB Fig. 1 fur die Bestimmung von MeB-Laufzeitdifferenzen erweitertes 
Blockschaltbild. 

[0031] Bei dem nachfolgend anhand des Blockschaltbilds in Fig. 1 beschriebenen Verfahren zum Bestimmen der 
Zielposition eines Schall ins Wasserabstrahlenden Ziels von einemzielfernen, untergetauchten Schallempfanger aus, 
wird dieTatsacheausgenutzt, daB sich im Wasser aufgrundvonTemperaturschichtungen oderaufgrund nurbegrenzter 
Wassertiefe Schallkanale bilden, in denen sich der Schall auf verschiedenen, reflexionsbedingten Ausbreitungswegen 
unterschiedlicher Lange ausbreitet und damit nach unterschiedlichen Laufzeiten am Schallempfanger eintrifft, so daB 
ein vom Schallempfanger aufgenommenes Empfangssignal zeitliche Uberlagerung des gleichen Signals enthalt. Ein 
Ziel ist jede im Wasser befindliche Schallquelle, wie sie z.B. von den Antriebsmaschinen oder sonstigen gerauscher- 
zeugenden Arbeitsaggregaten eines Oberflachenschiffes, U-Boots oderTorpedos gebildet wird. Der von diesen Ma- 
schinen und Aggregaten ausgesendete Schall ist breitbandig. Die Schallquelle kann aber auch ein aktiv sendendes 
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Sonar sein, das in einem Oberflachenschiff, U-Boot Oder Torpedo installiert 1st und Schalllmpulse ins Wasser abstrahlt. 
Der Schallempfanger 1 0, der auf einer ins Wasser eintauchenden Plattform, z.B. einem U-Boot, installiert ist, weist in 
bekannter Weise eine Mehrzahl von elektroakustischen Wandlern 11 auf, die im Ausfuhrungsbeispiel der Fig. 1 eine 
Linearantenne biiden, wie sie z.B. als Schleppantenne (Towed-Array) Oder als am Bootskorper befestigte Seitenan- 
tenne (Flank-Array), mitunter auch als Bordwandstreamer bezeichnet, bei U-Booten verwendet wird. 
[0032] Die elektrischen Ausgangssignale der Wandler 11 werden in bekannter Weise mittels eines Beamformers 
oder Richtungsbildners 12 richtungsselektiv zu Gruppensignalen addiert. Hierzu werden im Richtungsbildner 12 die 
Ausgangssignale der Wandler 11 zeit- oder phasenverzogert, und zwar derart, da(3 unter Berucksichtigung der ge- 
wiinschten Einf alls- oder Peilrichtung 8 m alle Empfangssignale der Wandler 1 1 konphas sind. Die Verzogerungszeiten 
x n m werden f ur jeden Wandler 1 1 fur eine vorgegeben Peilrichtung 6 m von dem Verzogerungsrechner 23 zur Verfiigung 
gestellt. Die in jeder Peilrichtung 6 m erhaftenen konphasen Ausgangssignale werden im Richtungsbildner 12 addiert 
und in einen Speicher 13 in Zuordnung zu der jeweiligen Peilrichtung 6 m abgelegt. Ein Pegeldetektor 14 ermittelt aus 
den Pegeln der abgespeicherten Gruppensignale die grbBten Pegel und gibt die den groBten Pegeln der Gruppensi- 
gnale zugeordneten Peilrichtungen 6 m als Zielpeilungen 8 Z aus, die einer Anzeigevorrichtung 15 zugefuhrt und in 
dieser numerisch und grafisch dargestellt werden. 

[0033] Fur eine Zielpeilung mit dem Peitwinkel 6 Z wird aus dem Speicher 1 3 das zugeordnete Gruppensignal aus- 
gelesen, und zur Bestimmung derZielposition des unter dem Peitwinkel &z detektierten Ziels wird dieses Gruppensignal 
einer Signalverarbeitung wie folgt unterzogen: 

[0034] In dem Gruppensignal werden mittels einer Laufzeitanalyse in dem Signalverarbeitungsblock 1 6 die im Grup- 
pensignal enthaltenen Laufzeitdifferenzen zwischen dem direkten Ausbreitungsweg des Schalls und den durch Refle- 
xionen im Schallkanal entstehenden indirekten Ausbreitungswegen des Schalls zwischen dem Ziel und dem Schall- 
empfanger 1 0 als MeB-Laufzeitdiff erenzen \q k gemessen . I n einem Modellblock 1 7 werden mittels eines ausgewahlten 
Modells der Schallausbreitung in einem Schallkanal fur eine angenommene Zielposition, die bestimmt ist durch eine 
angenommene Zielentfernung x und eine angenommene Zieltiefe z s , die Laufzeitdifferenzen zwischen dem direkten 
Ausbreitungsweg und ausgewahlten indirekten Ausbreitungswegen des Schalls vom Ziel zum Schallempfanger 1 0 als 
Modell-Laufzeitdifferenzen ty ( bestimmt. Die MeB-Laufzeitdifferenzen aus dem Signalverarbeitungsblock 16 und die 
Modell-Laufzeitdifferenzen aus dem Modellblock 17 werden einer Vergleichslogik 18 zugefuhrt. In der Vergleichslogik 

1 8 werden die Modell-Laufzeitdifferenzen mit dem MeB-Laufzeitdifferenzen verglichen. Ergibt der Vergleich einen nicht 
ausreichenden Grad der Ubereinstimmung, so wird die Bestimmung der Modell-Laufzeitdifferenzen im Modellblock 1 7 
mit einer geanderten Zielpositionsannahme wiederholt, und die neuen Modell-Laufzeitdifferenzen werden emeut mit 
den MeB-Laufzeitdifferenzen verglichen. Hierzu wird jeweils ein Zielpositions-Vorgabemodell 22 von dem Komparator 

19 aktiviert, der die neue Zielpositionsannahme der Vergleichslogik 18 und einem Tor 21 zugefOhrt. Dieser Vorgang 
wird soiange wiederhott, bis der Vergleich der Modell-Laufzeitdifferenzen mit den MeB-Laufzeitdifferenzen einen aus- 
reichenden Grad der Ubereinstimmung ergibt. Ist dies der Fall, dann wird von dem Komparator 1 9 das Tor 21 geoffnet 
und die zum Erreichen dieses Ubereinstimmungsgrads zuletzt getroffene Zielpositionsannahme als echte Zielposition 
mit der Zielentfernung x z * und der Zieltiefe z s * an die Anzeigevorrichtung 1 5 ausgegeben und in der Anzeigevorrichtung 
15 grafisch und numerisch dargestellt. 

[0035] Die Vergleichslogik 18, die durch Vergleich der Modell- und MeB-Laufzeitdifferenzen die Aufenthaltswahr- 
scheinlichkeit des Ziels an der angenommenen Zielposition ermittelt, kann nach verschiedenen Algorfthmen arbeiten. 
Im Ausfuhrungsbeispiel arbeitet die Vergleichslogik 1 8 nach dem Least-Square-Verf ahren . Die Wahrscheinlichkeit, bei 
der angenommenen Zielposition das Ziel zu finden, ist danach umgekehrt proportional zur Summe der kleinsten Ab- 
standsquadrate der Modell-Laufzeitdifferenzen und der MeB-Laufzeitdifferenzen. Dabei wird im Modell zunachst die 
Anzahl n M der zu berucksichtigenden Laufzeitdifferenzen festgelegt, die vom Signal-Noise-Ratio (SNR) des Gruppen- 
signals, von der maximalen Anzahl der indirekten Ausbreitungswege u.a. abhangt. Dann wird fur jede Modell-Lauf- 
zeitdifferenz t Mi diejenige Laufzeit t^i) gesucht, die einer beliebigen MeB-Laufzeitdiff erenz tc k am nachsten kommt. 
i ist der Index der Modell-Laufzeitdifferenzen und k der Index der MeB-Laufzeitdifferenzen. Der mathematische Aus- 
druck hierfur ist 

W0=t k (4) 

mit 

k = argmin ItMj-tcJ (5). 
[0036] Die Funktion argmin liefert das Argument, d.h. den Index k fiir den der Ausdruck auf der rechten Seite von 
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Gl.(5) minimal ist, hier also den Index k einer MeB-Laufzeltdifferenz mit minimalem Abstand zu einer Modell-Laufzeit- 
differenz. 

[0037] SchlieBlich wird die Summe der minimalen Abstandsquadrate uber alle Modell-Laufzeitdifferenzen gemaB: 

K(z s ,x) = f](t M(i - tjilf (6) 
gebildet. Das Minimum von K(z s x) liefert Zielentfemung und Zieltiefe. 

[0038] Die Summe K(z s ,x) dieser kleinsten Abstandsquadrate wird in einem Komparator 1 9 mit einem Konfidenzin- 
tervall Kl verglichen. Liegt K(z Sx ) nicht innerhalb dieses Konfidenzintervalls Kl, so wird vom Komparator 1 9 im Block 
20 eine neue Zielpositionsannahme mit einer neuen Zieltiefe z s und Zielentfemung x ausgelost, die dem Modellblock 
17 zugefuhrt wird und zu einer neuen Berechnung der Modell-Laufzeitdifferenzen fur die neue Zielpositionsannahme 
fuhrt. In der Vergleichslogik 18 wird wiederum die Summe K(z s x) der kleinsten Abstandsquadrate gebildet und im 
Komparator 19 mit dem Konfidenzintervall Kl verglichen. Liegt schlieBlich bei einer gewahlten Zielpositionsannahme 
die Summe K(z s x) der kleinsten Abstandsquadrate innerhalb des Konfidenzintervalls Kl, so wird von dem Komparator 
1 9 das Tor 21 geoff net, und die zuletzt getroffene Zielpositionsannahme wird als echte Zielposition mit der Zielentfer- 
nung x z * und Zieltiefe z s * der Anzeigevorrichtung 15 zugefuhrt und dort grafisch und numerisch dargestellt. 
[0039] Alternativ zu der Abprufung mit einem Konfidenzintervall konnen die Summen der minimalen Abstandsqua- 
drate in Wahrscheinlichkeiten umgerechnet werden und ein 2D- oder 3D-Plot fur die Aufenthaltswahrscheinlichkeit 
bestimmt wird. Ein MaB fur die Aufenthaltswahrscheinlichkeit P(z s x) eines Ziels in der angenommenen Zielposition 
ist der Kehrwert der Summe K(z s x) der minimalen Abstandsquadrate. Je geringer K(z s x) ist, desto groBer ist die 
Wahrscheinlichkeit, daB sich das Ziel bei den angenommenen Koordinaten x,z s der Zielposition befindet. Die Summe 
muB noch derart normiert werden, daB das raumliche Integral uber alle Wahrscheinlichkeitsdichten bzw. die Summe 
iiber alle P(z s x) gleich 1 ist, also gilt 



J Jp(z a ,x)dz, dx m 1 (7) 



[0040] Nach dieser Normierung konnen fur die Vielzahl der angenommenen Zielpositionen die Wahrscheinlichkeits- 
dichten pro Flache im Zielgebiet dargestellt werden. Ein Algorithmus oder ein Operator kann in der Darstellung erken- 
nen, ob eine mehrdeutige Situation mit mehreren lokalen Maxima und potentiellen Zielpositionen vorliegt oder ob die 
Situation eindeutig ist. 

[0041] Urn den EinfluB von stark gedampften, durch die Laufzeitanalyse nicht eindeutig identifizierbare MeB-Lauf- 
zeitdifferenzen hoherer Ordnung zu reduzieren, werden in der Vergleichslogik die minimalen Abstandsquadrate vor 
Summenbildung mit einem Wichtungsfaktor w, gewichtet, der in dem Block 22 generiert wird. Die Wichtungsfaktoren, 
die abhangig sind von der Anzahl der Schallreflexionen in indirekten Ausbreitungswegen, konnen direkt aus dem Schal- 
lausbreitungsmodell abgelegt werden. Bei einem nachfolgend noch beschriebenen einfachen Schallausbreitungsmo- 
dell finden beispielsweise Reflexionen des Schalls mit dem Reflexionsfaktor R in jedem indirekten Ausbreitungsweg 
statt, wobei die Amplitude der Schallwelle sich urn den Faktor R n reduziert. Da der Reflextionsfaktor R B am Boden und 
der Reflexionsfaktor R Q an der Oberflache des Schallkanals unterschiedlich sind, setzen sich die Wichtungsfaktoren 
W; in jedem Ausbreitungsweg i aus dem Teil mit Bodenreflexion und aus dem Teil der Oberflachenreflexion gemaB 

w i = R B B - R o n ° W 

zusammen, wobei n B die Anzahl der Bodenreflexionen und n 0 die Anzahl der Oberflachen reflexionen ist. Fur die Re- 
flexionsfaktoren R B und Rq werden plausible Werte zwischen 0 und 1 verwendet, die a priori festgelegt werden. Bei- 
spielsweise wird der Reflexionsfaktor R B = 0,5 und der Reflexionsfaktor R D = 0,75 festgelegt. Bei zwei Bodenreflexio- 
nen und einer Oberflachenreflexion im q-ten Ausbreitungsweg erga.be sich ein Wichtungsfaktor w q = 0,1875 mit dem 
das minimale Abstandsquadrat (ty q - t^fa)) 2 des Ausbreitungswegs i=q multipliziert werden muBte. 
[0042] Die Laufzeitanalyse zur Ermittlung der im Gruppensignal enthaltenden MeB-Laufzeitdifferenzen kann auf ver- 
schiedene Weise vorgenommen werden. In einem ersten Ausfuhrungsbeispiel wird zur Laufzeitanalyse die Autokor- 
relationsfunktion des Gruppensignals gebildet. Eine solche Bildung der Autokorrelationsfunktion ist in [2] beschrieben. 
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Die Autokorrelationsfunktion besteht bei Mehrwegeausbreitung aus einer Uberlagerung von Autokorrelationsfunktio- 
nen, die jeweils um die Differenzen der relativen Laufzeiten verschoben sind. Diese Laufzeitdifferenzen konnen durch 
Auswerten der Maxima der Autokorrelationsfunktion des Gruppensignals bestimmt werden, wobei aufgrund der Sym- 
metrie der Autokorrelationsfunktion nur positive Laufzeitdifferenzen berucksichtigt werden. In Fig. 2 ist eine solche 

s Autokorrelationsfunktion dargestellt. Deutlich ist ein Maximum bei einer Laufzeitdifferenz 0 ms zu sehen, das dem 
direkten Schallausbreitungsweg entspricht. Weitere Maxima sind bei 2 ms, 8 ms und 1 0 ms zu finden, die Laufzeitdif- 
ferenzen zwischen dem direkten Ausbreitungsweg und zwei indirekten Ausbreitungswegen sowie zwischen den beiden 
indirekten Ausbreitungswegen definieren. Die Maxima sind mit At 1 2 , At, 3, und At 23 gekennzeichnet. 
Bevorzugt wird zur Laufzeitanalyse im Signalverarbeitungsblock 1 6 eine Cepstralanalyse auf das Gruppensignal an- 

10 gewendet und die MeB-Laufzeitdifferenzen durch Auswertung der Maxima im Cepstrum gewonnen. Die Cepstralana- 
lyse ist an sich bekannt und in [3] beschrieben. Wie bereits erwahnt ist das reelle Cepstrum definiert als Fouriertrans- 
formierte des Logarithmus des Absolutwerts des Autospektrums des Gruppensignals gemaB Gl.(2), wobei X(e kot ) die 
Fouriertransformierte des Signals darstellt. Die Laufzeitdifferenzen lassen sich in den Maxima des Cepstrums leicht 
identifizieren. Ein Beispiel eines aus einer Cepstralanalyse gewonnenen Cepstrums ist in Fig. 3 dargestellt. Die Lauf- 

'5 zeitdifferenz At, ^ zwischen dem direkten Ausbreitungsweg und einem ersten indirekten Ausbreitungsweg ist bei 2 ms 
zu finden. Die Laufzeitdifferenz At 13 zwischen dem direkten Ausbreitungsweg und einen zweiten indirekten Ausbrei- 
tungsweg betragt 1 0 ms, wahrend noch eine Laufzeitdifferenz At 2 3 zwischen dem ersten und zweiten indirekten Aus- 
breitungsweg durch das Maximum bei 8 ms identifiziert wird. Ein Vergleich mit der Autokorrelationsfunktion gemaB 
Fig. 2 zeigt, daB diese Laufzeitdifferenzen in gleicher GroBe auch in der Autokorrelationsfunktion zu finden sind. Im 

20 Cepstrum lassen sich die Maxima aber wesentlich leichter identifizieren, weil hier keine Nebenmaxima auftauchen - 
wie dies in der Autokorrelationsfunktion der Fall ist - und somit reale Maxima von vorhandenen Reflexionen nicht mit 
Nebenmaxima verwechselt werden konnen. 

[0043] Als Modelle fur die Erstellung der Modell-Laufzeitdifferenzen konnen von dem bereits bekannten Modellen 
entsprechend geeignete ausgewahtt und in dem vorgestellten Verfahren eingesetzt werden. In [1] oder [5] sind ver- 

25 schiedene Modelle angegeben, die ggf. mit entsprechender Anpassung oder Modifikation verwendet werden konnen. 
Ein ebenfalls geeignetes Schallausbreitungsmodell ist unter dem Namen "RAY" bekannt und in [4] beschrieben. Aber 
auch die Verwendung eines einfachen geometrischen Schaliausbreitungsmodells, das im folgenden beschrieben wird, 
zeigt bereits recht gute Ergebnisse in der Genauigkeit der Bestimmung einer Zielposition. 
[0044] In Fig. 4 ist ein in einem Flachwassergebiet von z.B. einer Wassertiefe von ca. 200 m sich bildender Schall- 

30 kanal, in dem ein Ziel S in einer Zieltiefe z s und ein Schallempfanger E in der Empfangstiefe z E angeordnet sind, und 
voneinander eine Entfemung x aufweisen, schematisch dargestellt. Mit O ist die Wasseroberf lache und mit B der Boden 
des Schallkanals bezeichnet. Der Boden B wird als eben und das Schallgeschwindigkeitsprofil des Schallkanals als 
konstant mit der konstanten Schallgeschwindigkeit c angenommen. Der vom Ziel S abgestrahlte Schall gelangt einmal 
auf direktem Ausbreitungsweg zum Empfanger E und zum anderen auf indirekten Ausbreitungswegen durch Refle- 

35 xionen an der Oberflache O und am Boden B zu dem Empfanger E. Der direkte Ausbreitungsweg ist in Fig. 4 mit der 
Ordnungszahl k = 1 bezeichnet und ausgezogen dargestellt. Die indirekten Ausbreitungswege sind in Fig. 4 mit k = 2 
bis k= 5 bezeichnet, wobei der strichlinierte Ausbreitungsweg k= 2 eine Reflexion des Schalls an der Wasseroberf lache 
O und der strichliniert dargestellte indirekte Ausbreitungsweg k = 3 eine Reflexion des Schalls am Boden B enthalt. 
Die beiden weiteren strichpunktiert und punktiert dargestellten indirekten Ausbreitungswege k = 4 und k = 5 enthalten 

40 jeweils eine Reflexion des Schalls an der Wasseroberflache O und am Boden B. 

[0045] Der Weg, den der Schall in jedem der k-Ausbreitungswege zurucklegen muB und t k • c betragt, kann durch 
Annahme entsprechender Spiegelziele S' und S" oder Spiegelempfanger E' und E" fur die verschiedenen Ordnungen 
der Ausbreitungswege berechnet werden. Die Geometrie fur die Berechnung des indirekten Ausbreitungsweges k = 
4 ist in Fig. 4 durch das dargestellte rechtwinklige Dreieck, dessen Hypothenuse sich zwischen dem Spiegelziel S' und 

45 dem Spiegelempfanger E" erstreckt und dessen eine Kathete von der Entfemung x gebildet wird, herausgehoben. 
Allgemein berechnet sich die Laufzeit \ fur jeden der k-Ausbreitungswege gemaB Gl. (2) zu 



wobei z k der Weg ist, den der Schall ausschlieBlich in vertikaler Richtung zurucklegen muB. In Fig. 4 ist die Berechnung 
des Wertes z k fur die Ausbreitungwege k = 1 bis k = 5 angegeben. Der Elevationswinkel <t> k der Schalleinfallsrichtung 
am Empfanger E bei Schallausbreitung uber die Ausbreitungswege k berechnet sich zu 




(3), 



55 



(8) 
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und kann zur Korrektur des gemessenen Peilwinkels 8 Z verwendet werden. Der Einfallswinkel des uber den Ausbrei- 
tungsweg k = 4 am Empfanger E eintreffenden Schalls ist in Fig. 4 mit G> 4 eingezeichnet. 

Mittels des geometrischen Modells werden mit Vorgabe einerangenommenen Zielentfemung x und einerangenommen 
Zieltiefe z s die Laufzeiten t k fur ausgewahlte Ausbreitungswege des Schalls im Schallkanal bestimmt. Die Laufzeitdif- 

s ferenzen t^, , werden durch Subtraktion der Laufzeiten 1, fiir den direkten Ausbreitungsweg k = 1 von dem jeweiligen 
indirekten Ausbreitungsweg ^ mit k groBer 1 gewonnen und - wie beschrieben - der Vergleichslogik 1 8 zugefuhrt. 
[0046] In dem Fall, in dem zur Gewinnung der MeB-Laufzeitdifferenzen im Signalverarbeitungsblock 1 6 mittels Cep- 
stralanalyse das Cepstrum des Gruppensignals gebildet wird, kann das im Modellblock 1 7 verwendete Modell dahin- 
gehend erweitert werden, daft mit den Modell-Laufzeitdifferenzen ein Modell-Cepstrum gebildet wird. Die Vergleichs- 

'° logik 1 8 wird dahin abgewandelt, daft nun nicht mehr die einzelnen MeB-Laufzeitdifferenzen von Cepstrum und Modell 
miteinander verglichen werden, sondern unmittelbar das gemessene Cepstrum mit dem Modell-Cepstrum. 
[0047] Bei Inkaufnahme eines erweiterten Rechenaufwands kann auf die gesonderte Zieldetektion und Zielpeilung 
mittels der zur Fig. 1 beschriebenen Peilanlage verzichtet werden und mit dem erfindungsgemaBen Verfahren auch 
die Zieldetektion und Zielpeilung durchgefuhrt werden. Eine gegenuber Fig. 1 modifizierte Schaltungsanordnung fiir 

'5 das so abgewandelte Verfahren ist im Blockschaltbild in Fig. 5 dargestellt. In der Schattungsanordnung sind Gruppen- 
signal-Speicher und Pegeldetektor entfallen und der vorgegebene Wlnkel 8 m fiir den momentanen Schalleinfall (Peil- 
richtung) ist zusammen mit den Werten fur Zielentfemung x und Zieltiefe z s der Zielpositionsannahme an das Tor 21 
gelegt. Im iibrigen lauft das Verfahren in gleicher Weise ab, wie es zur Fig. 1 beschrieben ist, wobei durch Anderung 
der Peilrichtung 8 m in diskreten Winkelschritten der gesamte Azimut nach Zielen abgesucht wird. Sobald das ge- 

20 wunschte MaB der Ubereinstimmung zwischen den MeB-Laufzeitdifferenzen und den Modell-Laufzeitdifferenzen bei 
einer Peilrichtung 8,,, erreicht wird und der Komparator 1 9 das Tor 21 offnet, wird auch die zu der echten Zielposition 
gehorende Peilrichtung 8 m , die dann dem echten Peilwinkei 8 Z entspricht, der Anzeigevorrichtung 15 zugefuhrt und 
dort zusammen mit der zur echten Zielposition gehorenden echten Zielentfemung x* und echten Zieltiefe z s * numerisch 
und grafisch dargestellt. 

25 [0048] Das vorstehend beschriebene Verfahren laBt sich auch zum Bestimmen weiterer Zielparameter heranziehen, 
nachdem die Zielposition festgestellt worden ist. Diese Zielparameter sind die Frequenz des vom Ziel abgestrahlten 
Schalls, im folgenden Zielfrequenz genannt, und die radiale Geschwindigkeit des Ziels, das ist die Geschwindigkeits- 
komponente, mit dem sich das Ziel auf direktem, d.h. radialem, Weg auf den Empfanger zubewegt, im folgenden 
radiale Zielgeschwindigkeit genannt. Diese zusatzlichen Verfahrensschritte werden anhand des in Fig. 6 wiedergege- 

30 benen Blockschaltbilds erlautert, wobei das Blockschaltbild gem&8 Fig. 1 vollstandig in Fig. 6 integriert ist und insoweit 
auch die gleichen Bezugszeichen aus Fig. 1 in Fig. 6 ubernommen worden sind. 

[0049] Das Blockschaltbild der Fig. 1 ist dabei zunachst so erweitert, daB die Laufzeiten t k des Schalls, die mit dem 
Schallausbreitungsmodell im Modellblock 17 fur die Ausbreitungswege berechnet werden, fortlaufend in einen Lauf- 
zeitenspeicher 24 eingeschrieben werden, wobei die fur eine vorhergehende Zielposition berechneten Laufzeiten t^ 

35 jeweils mit den fur eine nachfolgende Zielposition berechneten Laufzeiten \ iiberschrieben werden, so daB im Lauf- 
zeitenspeicher 24 immer die fur die zuletzt angenommene Zielposition berechneten Laufzeiten Xy. abgelegt sind. Ist die 
wahre Zielposition erkannt worden, so wird durch das Ausgangssignal des Komparators 19 das Tor 21 geoffnet und 
die wahre Zielentfemung x z *, die wahre Zieltiefe z s * und der Peilwinkei B z der Anzeigevorrichtung 1 5 zugefuhrt. Gleich- 
zeitig gelangt das Komparatorsignal als Auslesesignal an den Laufzeitenspeicher 24, so daB die im Laufzeitenspeicher 

40 24 enthaltenen Laufzeiten ^ nacheinander ausgelesen und einem Fouriertransformator 25 als Steuersignal zugefuhrt 
werden. In dem zu Fig. 4 beschriebenen Ausfiihrungsbeispiel werden somit die Laufzeiten t k mit k=1 bis 5 ausgelesen. 
Im Fouriertransformator 25 wird das Gruppensignal, mit dem die Zielpeilung 8 Z festgestellt worden ist, einer Fourier- 
transformation unterzogen und furZeitpunkte, die den ausgelesenen Laufzeiten \ entsprechen, die Frequenzspektren 
gebildet. Aus den Frequenzspektren werden die Empfangsfrequenzen f Ek in Zuordnung zu den Laufzeiten t k ausge- 

45 lesen, und im Differenzbildner 26 werden die Differenzen der Empfangsfrequenzen 

AF k = f E1 -f Ek (9) 

50 mit k>1 gebildet und im Speicher 26 in Zuordnung zu den Laufzeiten \ mit k>1 abgespeichert. 

[0050] Die aus dem Laufzeitenspeicher 24 ausgelesenen Laufzeiten t* werden einem Rechner 28 zugefuhrt. Der 
Rechner 28 ist mit einem Vorgabemodul 29 fur die Zielfrequenz f 0 und radiale Zielgeschwindigkeit v R1 verbunden. Im 
Rechner 28 wird fiir jede Laufzeit t k die Dopplerverschiebung df k berechnet, wobei sich die Dopplerverschiebung fiir 
den direkten Ausbreitungsweg gemaB 
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c 



R1 



(10) 



und die Dopplerverschiebungen fur die indirekten Ausbreitungswege gemalJ 




(11), 



ergeben, wobei c die Schallgeschwindigkeit im Wasser und k die Ordnungszahl fur alle ausgewahlten indirekten Aus- 
breitungswege ist, die im Ausfuhrungsbeispiel der Fig. 4 gleich k=2, 3, 4, 5 betragt. Weiterhin werden im Rechner 28 
aus den berechneten Dopplerverschiebungen df k in den einzelnen Ausbreitungswegen kdie Differenzen der Doppler- 
verschiebungen gemaB 



bestimmt, wobei k wiederum die Ordnungszahl fur die indirekten Ausbreitungswege und damit k>1 ist. Diese Dopp- 
lerverschiebungs-Differenzen AF Mk werden in einem Speicher 30 in Zuordnung zu den Laufzeiten t k mit k>1 abgespei- 
chert. In einer Vergleichslogik 31 werden die Empfangsfrequenz-Differenzen AF k und die Dopplerverschiebungs-Dif- 
ferenzen AF Mk die den gleichen Schallaufzeiten t k mit k>1 zugehorig sind, miteinander verglichen. Ergibt der Vergleich 
einen nicht ausreichenden Grad der Ubereinstimmung, so wird die Bestimmung der Dopplerverschiebungs-Differenzen 
AF Mk im Rechner 28 mit einer geanderten Zielfrequenz f 0 und/oder einer geanderten radialen Zielgeschwindigkeit v R1 
wiederholt, und die neuen Dopplerverschiebungs-Differenzen AF^ werden wieder mit den Empfangsfrequenz-Diffe- 
renzen AF k verglichen. Dieser Vorgang wird solange wiederholt, bis der Vergleich der Dopplerverschiebungs-Differen- 
zen AF Mk mit den Empfangsfrequenzen-Differenzen AF einen ausreichenden Grad der Ubereinstimmung ergibt. Ist 
dies der Fall, dann wird ein Tor 32 geoffnet, an dem jeweils die zuletzt angenommene odervorgegebene Zielfrequenz 
f 0 und die zuletzt angenommene radiate Zielgeschwindigkeit v R1 anstehen. Mit Offnen des Tors 32 werden die zuletzt 
getroffenen Annahmen von Zielfrequenz f 0 und radialer Zielgeschwindigkeit v H1 der Anzeigevorrichtung 15 zugefuhrt 
und dort als echte Zielfrequenz f 0 * und echte radiale Zielgeschwindigkeit v R1 * graphisch und numerisch dargestellt. 
[0051] Die Vergleichslogik 31 arbeitet in gleicher Weise wie die zuvor beschriebenen Vergleichslogik 1 8, z,B. nach 
dem Least-Square-Verfahren, wobei als Kriterium beim Vergleich der Empfangsfrequenz-Differenzen AF k und der 
Dopplerverschiebungs-Differenzen AF Mk die Summe der minimalen Abstandsquadrate der Dopplerverschiebungs-Dif- 
ferenzen AF Mk zu ihren nachsten Nachbam in der Menge der Empfangsfrequenz-Differenzen AF k verwendet wird. Wie 
der Vergleichslogik 1 8 kann auch der Vergleichslogik 31 noch ein hier nicht dargestellter Komparator nachgeordnet 
sein, in dem die Summe derminimalen Abstandsquadrate mit einem Konfidenzintervall verglichen wird. Liegt die Sum- 
me innerhalb des Konfidenzintervalls, so ist der vorgegebene Grad der Ubereinstimmung der Dopplerverschiebungs- 
Differenzen AF Mk mit den Empfangsfrequenzen-Differenzen AF k erreicht, und der Komparator generiert das Toroff- 
nungssignal fur das Tor 32. Liegt die Summe auBerhalb des Konfidenzintervalls, so lost das Komparatorsignal eine 
neue Vorgabe von Zielfrequenz und/oder radialer Zielgeschwindigkeit im Vorgabemodul 29 aus. 
[0052] Urn eine groBere Anzahl von Vergleichswerten zu erhalten, ohne die Zahl der ausgewahlten Ausbreitungs- 
wege bzw. der ausgewahlten Laufzeiten t k in den Ausbreitungswegen zu vergroBern, werden nicht nur die Differenzen 
der Empfangsfrequenzen zwischen den Empfangsfrequenzen mit k>1 , die den ausgewahlten Laufzeiten t^ in den 
indirekten Ausbreitungswegen zugeordnet sind, und der Empfangsfrequenz f E1 , die der Laufzeit t, des direkten Aus- 
breitungswegs zugeordnet ist, gebildet, sondern auch die Empfangsfrequenz-Differenzen zwischen den indirekten 
Ausbreitungswegen selbst gebildet, so z.B. die Empfangsfrequenz-Differenz zwischen den Empfangsfrenquenzen f E3 
und f E2 , die der Schallaufzeit t 3 uber den Ausbreitungsweg dritter Ordnung und der Laufzeit t2 im Ausbreitungsweg 
zweiter Ordnung zugehdrig sind. In entsprechender Weise werden auch die Dopplerverschiebungs-Differenzen AF Mk 
zwischen den indirekten Ausbreitungswegen gebildet und in der Vergleichslogik 31 mit den entsprechenden Emp- 
fangsfrequenz-Differenzen in der beschriebenen Weise verglichen. 

[0053] Fur eine Verbesserung der Bestimmung der Zielposition und der Ermittlung der Geschwindigkeit des Ziels 
mit Mustersignalen wird das Modell im Modellblock 17 in Fig. 1 fur die Bestimmung von Modell-Stauchungsfaktoren 
s k erweitert. Fur die Berechnung der Modell-Laufzeitdifferenzen t^ und der Modell-Stauchungsfaktoren s k wird Fig. 4 
hinzugezogen und beispielhafl der indirekte Ausbreitungsweg k=3 mit einer Reflexion am Boden betrachtet, nachdem 
das Ziel S nach einem ersten MeBzeitpunkt mit einer Geschwindigkeit v in Richtung auf den Empfanger E gefahren 
ist und zum nachsten MeBzeitpunkt die Position S* eingenommen hat. Die Laufzeit uber den direkten Ausbreitungsweg 



AF Mk = df k- df 1 



(12), 



EP 1 271 175 A1 



k=1 betragt: 



t^Ux-v.lf + z, 8 



(13) 



und fur den Ausbreitungsweg uber den Boden 



t 3 =y(x- v -t) 2 +z 3 



(14). 



[0054] Die Modell-Laufzeitdifferenz t M betragt 



t M = t 3 -t t = - [j(x - vt) 2 + 2l 2 - ^(x - vt) 2 + z 3 2 ] (15). 



[0055] Die zeitliche Anderung der Modell-Laufzeitdifferenz t^ nach der Zeit t wird im Modellblock 1 7 berechnet und 
liefert den Modell-Stauchungsfaktor s k fur den zweiten MeBzeitpunkt: 



■ A*- 

s>, t) = — ^ = - (x - vt: 

dt c 



- vt) 2 + Zj a ^(x - vt) 2 + z 3 * 



(16) . 



[0056] Die zum ersten MeBzeitpunkt ermittelte Position des Ziels x, z s wird eingesetzt und s k fur unterschiedliche 
Geschwindigkeiten v zusammen mit den Modell-Laufzeitdifferenz gemaB Gl. (15) t^, bestimmt. 
[0057] Die zeitliche Anderung der Modell-Laufzeitdifferenz gemaB Gl. (1 6) entspricht einer zeitlichen Anderung der 
Phase <p bei einer Frequenz <o=2jif 



d9 = (0 ^M 

dt at 



(17) 



[0058] Die zeitliche Anderung der Phase cp ruft eine Dehnung oder Stauchung des zeitlichen Verlaufs des Gruppen- 
signals hervor, die durch die Geschwindigkeit des Ziels entsprechend dem Dopplereffekt verursacht wird, und fuhrt zu 
einer Dopplerverschiebung Aw und einer Drehung der Phase <p des komplexen Frequenzspektrums bei jeder Frequenz 
w Oder des zeitlichen Verlaufs des empfangenen Gruppensignals langs jedem Ausbreitungsweg. Die Dopplerverschie- 
bung Aoabezogen auf die Frequenz to istgleich dem Modell-Stauchungsfaktor s k und abhangig von der Geschwindigkeit 
v des Ziels sowie der Zeit t, in der sich das Ziel mit der Geschwindigkeit v bewegt. 



(18) 



[0059] Ein Gruppensignal, das einen Impuls enthalt, erhalt uber die indirekten Ausbreitungswege k>1 ebenfalls eine 
relative Stauchung Oder Dehnung der ausgesendeten Impulslange, die durch den Modell-Stauchungsfaktor s k beruck- 
sichtigt ist. Mit dem Modell-Stauchungsfaktoren s k wird das zum ersten MeBzeitpunkt empfangene Zeitstuck des Grup- 
pensignals gestaucht und/oder gedehnt und dadurch Mustersignale hergestellt. 

[0060] Fig. 7 zeigt ein gegenuber Fig. 1 erweitertes Blockschaltbild. Der Signalverarbeitungsblock 16 fur die Lauf- 
zeitanalyse ist wahlweise durch Offnen eines Schalters 161 uber einen BandpaB 162 der Mittenfrequenz w mit dem 
Speicher 13 fur die Gruppensignale am Ausgag des Richtungsbildners 12 verbunden. Das Gruppensignal wird nach 
Detektion eines Ziels einer Autokorrelationsschaltung 163 im Signalverarbeitungsblock 16 zugefuhrt, in der die Auto- 
korrelationsfunktion des Gruppensignals gebildet und Maxima bei den MeB-Laufzeitdifferenzen At 1 2 ; At^ 3 ; At 1 3 aus- 
gewertet werden, wie in Fig. 2 gezeigt. Die MeB-Laufzeitdifferenzen t^ werden in der Vergleichslogik 18 mit den 
Modell-Laufzeitdifferenzen t Mk verglichen und ergeben die Zielposition x, z s zum ersten MeBzeitpunkt. Das Gruppen- 
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signal wird auBerdem zum gieichen MeBzeitpunkt in einer Mustersignalschaltung 165 bezuglich seiner GroBe und 
Stabilitat untersucht und der zeitliche Verlauf des Gruppensignals mit dem groBten Pegel und der groBten Stabilitat 
innerhalb eines Zeitintervalls ATabgespeichert. In dem Modellbiock 17werden Muster der Modell-Stauchungsfaktoren 
s k uber den Modell-Laufzeitdifferenzen t M der zum nachsten MeBzeitpunkt angenommenen Zielpositionen und Ge- 

5 schwindigkeiten aus dem Zielpositions- und Geschwindigkeitsvorgabemodell 22 berechnet und bereitgestellt. Zur Er- 
mittlung von Mustersignalen wird das abgespeicherte Zeitstuck des Gruppensignals entsprechend den ermittelten 
Modell-Stauchungsfaktoren s k gestaucht oder gedehnt, beispielsweise durch unterschiedliche Auslesezeiten des ab- 
gespeicherten zeitlichen Verlaufs. Zum nachsten MeBzeitpunkt nach dem Zeitintervall AT wird das empfangene Grup- 
pensignal mit den Mustersignalen aus der Mustersignalschaltung 165 in einer Kreuzkorrelationsschaltung 166 korre- 

10 liert, dessen Integrationszeit-t die gleiche Lange wie die Integrationszeit in der Autokorrelationsschaltung 1 63 aufweist. 
Die bei der Kreuzkorrelation samtlicher Mustersignale mit dem Gruppensignal ermittelten Maxima werden zum Be- 
stimmen von MeB-Laufzeitdifferenzen t* ausgewertet und mit den Modell-Laufzeitdifferenzen ty in der Vergleichslogik 
18 verglichen. 

[0061] Bei Dbereinstimmung werden die in dem Zielpositions- und Geschwindigkeitsvorgabemodell 22 vorgegebe- 
'5 nen Werte fur die Zielposition x z *, z s * und die Geschwindigkeit v* uber das Tor 21 an die Anzeigevorrichtung 1 5 durch- 
geschaltet und angezeigt. 
[0062] Literatur: 

[1 ] Computational Ocean Acoustics, F. B. Jensen, W. A. Kuperman, M. B. Porter, H. Schmidt, AIP Press American 
20 Istitute of Physics, New York 1994, insbesondere S. 58, 

[2] J. S. Bendat and A. G. Piersol "Engineering Applications of Correlation an spektral Analysis", John Wiley and 
Sons, Inc. 1993. 

25 [3] B. Bogert, M. Healy and J. Tukey, 'The quefrency alanysis of time series for echoes: Cepstrum, pseudo-auto- 

covariance, cross-cepstrum and saphe cracking" Proc. Symposium Time Series Analysis, 1 963, pages 209 - 243. 

[4] Bowlin et al "Ocean Acoustics Ray-Tracing Software RAY", WHOI Technical Report, WHOI-93-10. (ftp;//oalib. 
saic.com/pub/oalib/ray) 

30 

[5] Sound Progation in the Sea, Robert J. Urick, Peninsula Publishing, 19B2, ISBN 0-9321 46-0B-2, insbesondere 
S. 8-11 

35 Patentanspruche 

1. Verfahren zum Bestimmen der Zielposition eines Schall ins Wasser abstrahlenden Ziels unter Auswertung des 
vom Ziel abgestrahlten Schalls, derauf direktem Ausbreitungsweg und auf durch Reflexionen entstehenden indi- 
rekten Ausbreitungswegen von einem untergetauchten, zielfernen Schallempfanger (1 0) mit einer Anzahl von elek- 

40 troakustischen Wandlern (1 1 ) empf angen wird, dadurch gekennzeichnet, daB die elektrischen Ausgangssignale 

derWandler(11) richtungselektiv zu Gruppensignalen addiert werden, daB in den Gruppensignalen mittels Lauf- 
zeitanalyse Laufzeitdifferenzen des Schalls zwischen dem direkten Ausbreitungsweg und indirekten Ausbreitungs- 
wegen als MeB-Laufzeitdifferenzen bestimmt werden, daB mittels eines Modells der Schallausbreitung fur eine 
angenommene Zielposition Laufzeitdifferenzen zwischen dem direkten Ausbreitungsweg und ausgewahlten indi- 

45 rekten Ausbreitungswegen vom Ziel zum Schallempfanger (10) bestimmt werden und diese Modell-Laufzeitdiffe- 

renzen mit den MeB-Laufzeitdifferenzen verglichen werden, daB die Bestimmung der Modell-Laufzeitdifferenzen 
mit jeweils geanderter Zielpositionsannahme solange wiederholt wird, bis der Verglebh der Modell-Laufzeitdiffe- 
renzen mit den MeB-Laufzeitdifferenzen einen vorgegebenen Grad der Ubereinstimmung ergibt, und daB die zum 
Erreichen dieses Obereinstimmungsgrads getroffene Zielpositionsannahme als echte Zielposition ausgegeben 

so wird. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, daB als Kriterium beim Vergleich von Modell- und 
MeB-Laufzeitdifferenzen die Summe der minimalen Abstandsquadrate der Modell-Laufzeitdifferenzen zu ihren 
nachsten Nachbarn in der Menge der MeB-Laufzeitdifferenzen verwendet wird. 
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3. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daB der vorgegebene Grad der Ubereinstimmung dann 
erreicht ist, wenn die Summe der minimalen Abstandsquadrate innerhalb eines Konfidenzintervalls liegt. 
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4. Verfahren nach Anspruch 2, dadurch gekennzeichnet, daB die Wahrscheinlichkeit des Aufenthalts des Ziels in 
der angenommenen Zielposition als Kehrwert der Summe der minimalen Abstandsquadrate ermittelt wird und da3 
aus der Vielzahl der fur angenommene Zielpositionen erhaltenen Wahrscheinlichkeiten die Wahrscheiniichkeits- 
dichte pro Flacheneinheit des Zielgebiets bestimmt wird und anhand der Wahrscheinlichkeitsdichte potentielle 
Zielpositionen aufgesuchtwerden. 

5. Verfahren nach einem der Anspruche 2 - 4, dadurch gekennzeichnet, daB die minimalen Abstandsquadrate vor 
der Summenbildung mit einem Wichtungsfaktor (w ( ) multipliziert werden. 

6. Verfahren nach Anspruch 5, dadurch gekennzeichnet, daB als Wichtungsfaktor (Wj) eine Potenz der Anzahl 
(n B n 0 ) der Reflexionen mit den Reflexionsfaktoren (R B , Rq) im Ausbreitungsweg (i) gemaB 



benutzt wird, wobei der Reflexionsfaktor (R B ) fur eine Reflexion am Boden und der Reflexionsfaktor (Rq) fur eine 
Reflexion an der Oberflache (O) eines durch Boden (B) und Oberflache (0) begrenzten Schallkanals plausibel 
zwischen 0 und 1 gewahlt wird. 

7. Verfahren nach einem der Anspruche 1 - 6, dadurch gekennzeichnet, daB zur Laufzeitanalyse die Autorkorre- 
lationsfunktion des Gruppensignals bestimmt wird und die MeB-Laufzeitdifferenzen durch Auswerten der Maxima 
der Autokorrelationsfunktion gewonnen werden. 

8. Verfahren nach einem der Anspruche 1 - 6, dadurch gekennzeichnet, daB zur Laufzeitanalyse eine Cepstral- 
analyse auf die Gruppensignale angewendet wird und daB die MeB-Laufzeitdifferenzen durch die Auswertung der 
Maxima im Cepstrum gewonnen werden. 

9. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, daB zur Laufzeitanalyse eine Cepstralanalyse auf die 
Gruppensignale angewendet und jeweils das Cepstrum des Gruppensignals als MeB-Cepstrum gebildet wird und 
daB mit den Modell-Laufzeitdifferenzen ein Modell-Cepstrum gebildet und der Vergleich der MeB-Laufzeitdifferen- 
zen und der Modell-Laufzeitdifferenzen durch unmittelbaren Vergleich des Modell-Cepstrums mit dem MeB-Cep- 
strum durchgefuhrt wird. 

10. Verfahren nach einem der Anspruche 1 - 9, dadurch gekennzeichnet, daB als Modell ein geometrisches Modell 
der Schallausbreitung in einem Flachwassergebiet mit der Wassertiefe (z w ) bei ebenem Gewasserboden und 
konstantem Schallgeschwindigkeitsprofil zugrunde gelegt wird, daB bei einer bekannten Tiefe (z E ) des Schallemp- 
fangers und einer angenommenen Zieltiefe (z s ) und einer angenommen azimutalen Zielentfernung (x) die Laufzeit 
(t k ) fur k ausgewahlte Ausbreitungswege gemaB 



berechnet wird, wobei c die konstante Schallgeschwindigkeit und z k der Weg ist, den der Schall in jedem der k 
Ausbreitungswege ausschlieBlich in vertikaler Richtung zurucklegt, und daB zur Gewinnung der Modell-Laufzeit- 
differenzen die Laufzeit (t^ fur den direkten Ausbreitungsweg von den Laufzeiten (y fur k ausgewahlte indirekte 
Ausbreitungswegen mit k> 1 subtrahiert wird. 

11. Verfahren nach einem der Anspruche 1-10, dadurch gekennzeichnet, daB die Pegel der Gruppensignale be- 
stimmt und durch Auswertung der groBten Pegel Ziele detektiert und die Zielpeilungen bestimmt werden und daB 
nur in denjenigen Gruppensignalen MeB-Laufzeitdifferenzen bestimmt und mit Modell-Laufzeitdifferenzen vergli- 
chen werden, in denen ein Ziel detektiert wurde. 

12. Verfahren nach einem der Anspruche 1-11, dadurch gekennzeichnet, daB mit Feststellung der echten Zielpo- 
sition mit dem Modell der Schallausbreitung bei der Bestimmung der Modell-Laufzeitdifferenzen berechnete Lauf- 
zeiten des Schalls fur den direkten Ausbreitungsweg und fur die ausgewahlten indirekten Ausbreitungswege zur 
Zielposition aufgerufen werden, daB aus dem Gruppensignal, in dem das Ziel detektiert wurde, Frequenzsprektren 
zu den den aufgerufenen Laufzeiten (t k ) entsprechenden Zeitpunkten gebildet werden, daB aus den Frequenz- 
spektren Empfangsf requenzen (f EK ) ausgelesen, die Differenzen der Empf angsfrequenzen (AF k ) gebildet und den 



n B i 
w, = R B • R 0 
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aufgerufenen Laufzeiten (y zugeordnet werden, daB fur eine angenommene, vom Ziel abgestrahlte Zielfrequenz 
(f 0 ) und fur eine angenommene radiale Zielgeschwindigkeit (v R ) zu jeder der aufgerufenen Laufzeiten (t k ) eine 
Dopplerverschiebung (df k ) berechnet und die Differenzen der Dopplerverschiebungen (df k ) gebildet werden, daB 
die Empfangsfrequenz-Differenzen (AF k ) und die Dopplerverschiebungs-Differenzen(AF Mk ), die den gleichen Lauf- 
zeiten (t k ) zugehorig sind, miteinander verglichen werden, daB die Berechnung der Dopplerverschiebungen (df k ) 
mit jeweils geanderter Zielfrequenz (f 0 ) und/oder radialer Zielgeschwindigkeit (v R ) und der Vergleich mit den Emp- 
fangsfrequenz-Differenzen (AF k ) solange wiederholt wird, bis sich ein vorgegebener Grad der Obereinstimmung 
ergibt, und daB die zum Erreichen dieses Ubereinstimmungsgrads getroffenen Annahmen von Zielfrequenz (f 0 ) 
und radialer Zielgeschwindigkeit (v R ) als wahre Zielfrequenz (f 0 *) und wahre radiale Zielgeschwindigkeit (v* R1 ) 
ausgeben werden. 

1 3. Verf ahren nach Anspruch 1 2, dadurch gekennzeichnet, daB als Kriterium beim Vergleich von Empfangsfrequenz- 
Differenzen (AF k ) und Dopplerverschiebungs-Differenzen (AF Mk ) die Summe der minimalen Abstandsquadrate 
der Dopplerverschiebungs-Differenzen (AF Mk ) zu ihren nachsten Nachbam in der Menge der Empfangsfrequenz- 
Differenzen (AF k ) verwendet wird. 

14. Verfahren nach Anspaich 13, dadurch gekennzeichnet, daB der vorgegebene Grad der Obereinstimmung dann 
erreicht ist, wenn die Summe der minimalen Abstandsquadrate innerhalb eines Konfidenzintervalls liegt. 

15. Verfahren nach einem der Anspruche 12 - 14, dadurch gekennzeichnet, daB die Dopplerverschiebung (df,) im 
direkten Ausbreitungsweg als Produktaus an genommener Zielfrequenz (f 0 ) und angenommener radialer Zielge- 
schwindigkeit (v R1 ) dividiert durch die Schallgeschwindigkeit (c) im Wasser und die Dopplerverschiebung (df k , k>1 ) 
fur die ausgewahlten indirekten Ausbreitungswege als die mit dem Verhaltnis von Laufzeit (t,) im direkten Aus- 
breitungsweg und Laufzeit^, k>1) im jeweiligen indirekten Ausbreitungsweg multiplizierte Dopplerverschiebung 
(df^ im direkten Ausbreitungsweg berechnet wird. 

16. Verfahren nach einem der Anspruche 12 -15, dadurch gekennzeichnet, daB als Empfangsfrequenzen (f Ek ) si- 
gnifikante Empfangsfrequenzbanderaus den Freuquenzspektren ausgelesen werden und die Zielfrequenz (f 0 ) als 
ein vom Ziel abgestrahltes Frequenzband angenommen wird. 

17. Verfahren nach Anspruch 1 6, dadurch gekennzeichnet, daB Harmonische der Zielfrequenz im Antriebsgerausch 
zur Klassifizierung des Ziels verwendet werden. 

18. Verfahren nach einem der Anspruche 1 - 11, dadurch gekennzeichnet, daB derzeitliche Verlauf des Gruppen- 
signals, in dem ein Ziel detektiert und eine Zielposition (x, z s ) zugeordnet ist, innerhalb eines Zeitintervalls (AT) 
gespeichert wird und zum Erstellen von Mustersignalen mit vorgebbaren Modell-Stauchungsfaktoren (s k ) gedehnt 
und/oder gestaucht wird, daB derModell-Stauchungsfaktor(s k ) fur jede Modell-Laufzeitdifferenz der zum nachsten 
MeBzeitpunkt angenommenen Position aus der ersten Zielposition und vorgebbaren Werten fur die Geschwindig- 
keit des Ziels bestimmt wird, daB Modell-Stauchungsfaktor (S k ) den Modell-Laufzeitdifferenzen (t Mi ) zugeordnet 
werden, daB die Mustersignale mit dem zum nachsten MeBzeitpunkt in einem Zeitintervall empfangenen Grup- 
pensignal kreuzkorreliert werden und daB die MeB-Lauf zeitdifferenzen (t k ) durch Auswerten der Maxima der Kreuz- 
korrelationsfunktion gewonnen werden, mit den Modell-Laufzeitdifferenzen (t^) verglichen werden und Position 
und Geschwindigkeit des Ziels liefem. 

19. Verfahren nach Anspruch 18, dadurch gekennzeichnet, daB die erste Zielposition mit den Me(3-Laufzeitd'rfferen- 
zen (y aus der Auswertung der Autokorrelationsfunktion ermittelt wird. 

20. Verfahren nach Anspruch 18 oder 19, dadurch gekennzeichnet, daB eine zeitabhangige relative Dopplerver- 
schiebung langs der Ausbreitungswege aus den zu den ermittelten MeB-Laufzeitdifferenzen gehorenden Modell- 
Stauchungsfaktoren (s k ) bestimmt wird. 

21 . Verfahren nach Anspruch 20, dadurch gekennzeichnet, daB aus der relativen Dopplerverschiebung (Aw/©) und 
der Geschwindigkeit v eine ZielWassifizierung hergeleitet wird. 

22. Verfahren nach Anspruch 21 , dadurch gekennzeichnet, daB Mittenfrequenz und Bandbreite des Gruppensignals 
ausgewah It werden. 

23. Verfahren nach Anspruch 22, dadurch gekennzeichnet, daB die Auswahl an die vom Ziel zu erwartende Schall- 
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abstrahlung angepaBt wird. 

24. Verfahren nach Anspruch 20 oder21 , dadurch gekennzeichnet, daB die Anpassung abhangig von der Entfernung 
zum Ziel und/oder zeitabhangig gewahlt wird. 

25. Verfahren nach einem der Anspruche 23 - 24, dadurch gekennzeichnet, daB die Bandbreite an eine gewunschte 
Entfernungsauflosung angepaBt wird. 

26. Verfahren nach einem der Anspruche 1 8 - 25, dadurch gekennzeichnet, daB die Integrationszeit (i) fur die Kreuz- 
korrelationsfunktion und die Geschwindigkeitsauflosung aufeinander abgestimmt werden. 
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VC^X 2 ^ 2 

fur k = 1: Zt =Z S -Z E (Direkt) 

k = 2: Z2 = Z S + Z£ (0) 

k = 3: Z 3 = 2Z W -Z S -Z E (B) 

k = A: Z A = 2Z W + Z S -Z E (OB) 

k = 5: Z 5 = 2Z W -Z S + Z E (BO) 

Fig. 4 
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